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Abstract 

Motivated by the conjectured Penrose inequality and by the work of Hawking, 
Geroch, Huisken and Ilmanen in the null and the Riemannian case, we examine 
necessary conditions on flows of two-surfaces in spacetime under which the Hawk- 
ing quasilocal mass is monotone. We focus on a subclass of such flows which we 
call uniformly expanding, which can be considered for null as well as for space- 
like directions. In the null case, local existence of the flow is guaranteed. In the 
spacelike case, the uniformly expanding condition leaves a 1-parameter freedom, 
but for the whole family, the embedding functions satisfy a forward-backward 
parabolic system for which local existence does not hold in general. Neverthe- 
less, we have obtained a generalization of the weak (distributional) formulation of 
this class of flows, generalizing the corresponding step of Huisken and Ilmanen's 
proof of the Riemannian Penrose inequality. 

1 Introduction 

Penrose has conjectured [1] that in an asymptotically flat spacetime satisfying the 
dominant energy condition, the total ADM mass M and the area \S\ of an outermost 
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apparent horizon S must satisfy 



S 



(1) 



M > 



16tt 



in units where Newton's gravitational constant is unity. A similar inequality should 
hold for the Bondi mass as well. In the presence of a negative cosmological constant 



is conjectured [2] for suitably defined masses at null infinity, where x is the Euler char- 
acteristic of S. Penrose gave a heuristic argument for (1) based on the standard view 
of gravitational collapse, involving cosmic censorship, a final stationary configuration 
and the classification of stationary black holes, among others. 

For inequality (1) to hold, it is known that S has to be area outer minimizing in 
some sense [3, 4]. For counterexamples when this condition is not fullfilled see [4, 5]. 
However, a rigorous proof of this Penrose inequality has been achieved only in two 
particular cases: in spherical symmetry [6, 7, 8] and in the time-symmetric case. The 
latter can be formulated as a purely Riemannian problem, termed the Riemannian 
Penrose inequality, which involves a Riemannian, asymptotically Euclidean 3-manifold 
(E, 7), with non-negative Ricci scalar, which translates the dominant energy condition, 
and with an outermost minimal surface S, which substitutes the apparent horizon. The 
inequality (1) for this case was proven by Huisken and Ilmanen [9] for connected S and 
by Bray [10], with a totally different method, for arbitrary S. Bray's method uses a 
flow of Riemannian metrics interpolating between the starting metric 7 and the metric 
of a t = const, slice of Schwarzschild outside the horizon. This flow of metrics keeps the 
area of the outermost minimal surface constant and does not increase the ADM mass. 
The Penrose inequality follows from the fact that equality holds for Schwarzschild. 

The approach of Huisken and Ilmanen makes use of the Geroch-Hawking mass [11] 
which is associated to any closed surface S in a hypersurface E: 



where rjs is the surface element, \S\ is the total area of S, and p is the mean curvature 
of S in E. The basic observations due to Geroch [11] and to Jang and Wald [12] refer 
to the case A = and are the following: (3) is monotonic for a special flow of two- 
surfaces called the Inverse Mean Curvature Flow (IMCF), approaches the ADM mass 
for suitable surfaces at infinity, such as metric 2-spheres, and equals the right hand 
side of (1) at the minimal surface. Thus, if the IMCF exists globally and approaches 
infinitely large round spheres, the Penrose inequality follows. However, the IMCF does 
in general develop singularities, and Huisken and Ilmanen [9] succeeded in defining a 
suitable weak version of the flow for which global existence could be proven. This weak 
version of the IMCF allows for jumps and one needs to question monotonicity of the 
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Geroch-Hawking mass at the jumps. The condition of connectedness of the starting 
minimal surface was required at this point. 

The Geroch-Hawking mass is a translation to 3-dimensional Riemannian manifolds 
of the Hawking mass [13] 

— * 

defined for 2 surfaces S embedded in spacetime (V, g) with mean curvature vector H, 
where H 2 = (H ■ H), the dot denoting the inner product on (V, g). US is embedded in 
a 3-surface E such that H points along the unique unit normal v of S, one can write 
H = pi?. However, a 2-surface S will in general have different Hawking and Geroch- 
Hawking masses, the former depending only on S and the latter also on E. We believe 
that the Hawking mass is a good candidate for proving the general Penrose inequality 
along the lines sketched above for Mq in the Riemannian case. In spherical symmetry, 
there is a quite simple and otherwise general proof based on monotonicity of Mh [14]. 
To summarize existing work, to describe the extensions achieved in the present paper 
and to formulate the open issues, we found it useful to distinguish the following five 
steps: 



1. Variation of the Hawking mass 

2. Monotonicity properties for special flows 

3. Local existence of good flows 

4. Global existence of good flows 

5. Limits at the horizon and at infinity. 

As to the variation of Mh, in Sect. 2 we find an expression which is valid for 
arbitrary variations and which, when the flow is everywhere spacelike, takes the form 



dM H {S x ) 



L [matter 



dX V I ' > . I Jg 



density + gravitational energy density + rest] r) Sx (5) 



where A denotes the flow parameter and S\ the corresponding two-surface. The names 
of the first two terms in the bracket mean that they vanish in vacuum and in spherical 
symmetry, respectively, and they turn out to be non-negative under reasonable positive- 
energy conditions on the matter, and under causality conditions on H and on the flow 
vector £. A "gravitational energy density" is not well-defined in general, but the point 
here is to perform some intelligent splitting which allows the "rest" (or its integral) to 
be removed by some simple and not over-restrictive choices of the flow. 

Such monotonicity conditions are specified in Theorem 1 in Section 3, which is the 
main result of this paper and which we sketch here. To have the energy terms non- 
negative, we impose the dominant energy condition, and we require that H is spacelike 
or null with (£ • H) > 0. To get rid of the "rest" we require two conditions on the 
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flow. The first one is the dual inverse mean curvature condition, reading that (£ • H*) 
should be a constant c(A) on each S\, where H* denotes the dual of H. For the second 
condition, our splitting suggests two natural alternatives. The first one is the inverse 
mean curvature condition, namely that (£ • H) — a(A) with a(A) constant on each S\. 
This condition guarantees, in particular, an IMCF on the three-surface spanned by the 
evolving two-surfaces themselves, so it is a natural extension of the Geroch condition 
mentioned above. To get monotonicity of Mh, we take the inverse mean curvature 
condition and the dual inverse mean curvature condition together, for which we adopt 
the name uniformly expanding flow. In this case, monotonicity still holds if we allow 
£ to be null somewhere, or everywhere. On regions where either £ or H are null, the 
inverse mean curvature condition in fact implies the dual condition, and vice versa. If 
in particular £ is null, this condition has been analyzed by Hawking himself [13], by 
Eardley [15] and, more extensively, by Hayward [6], while in the spacelike case, the 
uniformly expanding flows were investigated by Malec, Mars and Simon [16]. In the 
special case c(A) = 0, the flow is tangent to the inverse mean curvature vector and was 
analyzed by Frauendiener [17]. While a(A) can be set equal to 1 by a reparametrization, 
if it has no zeros, the presence of c(A) provides a one-parameter freedom for the flows 
in the spacelike case. 

The alternative to the inverse mean curvature condition arising from the present 
work is a Poisson equation on S which determines the lapse function of the flow, with 
a source term depending on the scalar curvature of S and on H 2 . This equation can 
always be solved uniquely on any S\. Here we restrict ourselves to flows which are 
spacelike everywhere. The condition seems new, and has a counterpart in codimension 
one, with H 2 replaced by p 2 , where it can be used as an alternative to IMCF in order to 
obtain monotonicity of the Geroch mass. The properties of this flow will be elaborated 
elsewhere. 

The conditions on the flow vector £ as mentioned above guarantee monotonicity 
and they are consistent and solvable on every two-surface S\. However, this does 
not mean that we have shown existence of a flow in spacetime, i.e. a sequence of S\, 
not even locally. For the IMCF in 3-dimensions, the embedding functions satisfy a 
parabolic system for which local existence is guaranteed. In spacetime, the uniformly 
expanding condition gives a local existence result in the case when £ is null, since only 
the null geodesies of the flow have to be determined in this case. On the other hand, 
in the spacelike case we are left with a forward-backward parabolic system which was 
observed by Huisken and Ilmanen in the case c(A) = and which we show in general 
in Sect. 4. For such systems, not even local existence results are available. 

A key element of the Huisken-Ilmanen proof in codimension 1 is, however, the weak 
or distributional formulation of the flow in terms of its level sets, which satisfy degen- 
erate elliptic equations, if sufficiently differentiable. We have obtained a generalization 
of this variational principle to the spacetime case, which we also describe in Sect. 4. 

We finally comment on point 5 of the list above: Mh obviously gives the right hand 
side of (2) on any marginally trapped surface. At infinity, provided the S\ approach 
surfaces of constant curvature, the Hawking mass approaches expressions which have 
been studied before: ADM and Bondi for A = 0; for A < see [2]. Except for the 
asymptotically flat case at spatial infinity, it is an open issue if the flow necessarily 
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leads to such surfaces. 



2 Variation of the Hawking mass 

We consider orientable Riemannian manifolds (V, g) and orientable closed (i.e. compact 
without boundary) submanifolds (S, h) , where h is the induced metric on S. For our 
main results, we have to restrict V to be 4-dimensional with g of signature (— , +, +, +), 
and S to be 2-dimensional, with a positive definite metric h. We will also consider the 
case that S is embedded in a 3-dimensional manifold £ with positive definite metric, 
which follows from the above one by choosing a suitable embedding £ C V. Some 
auxiliary results, in particular Lemma 1 and Lemma 2, hold for arbitrary dimensions 
and arbitrary metric signatures. Up to and including these Lemmas, we will consider 
this general setup, and impose suitable restrictions afterwards. 

Our aim in this section is to calculate the variation of the Hawking mass (4) along 
an arbitrary flow. We take an arbitrary C 1 embedding 

$ : S x / — ► V, 

(a;, A) — ► $(x,\) = $x(x), 

where S is a copy of S viewed as an abstract manifold detached from V and I 3 is an 
open interval of R. The leaves S\ = <&\(S) are required to be C 2 , so that the variation 
of H can be defined, and we write S for So. Choosing an orientation on (S, h) defines 
uniquely a canonical two-form t]s from the volume form r\ on (V,g). 

The tangent space T p of V at p G S can be decomposed as T p = T p © T p L where T p 
and T p L are the tangent and normal spaces to S, respectively, and we shall write v = 
tfll +V 1 - for any vector v E T p . The flow of submanifolds S\ has a flow vector v = d&(d\) 
which may have an arbitrary causal character. Accordingly, the embedded submanifold 
M = $(5 x /) need not have constant causal character at different points; in the 
spacetime case, the hypersurface Af could be timelike, spacelike or null at different 
points. Let us define £ as its normal component with respect to S\, £ = tr 1 . We shall 
call this vector the normal flow vector and also simply flow vector when no confusion 
arises. 

Evaluating the Hawking mass for any leaf S\, we write Mff(A) = Mh(S\). We will 
not add a subscript A to H or to other geometric objects on S\, unless the meaning is 
not clear from the context. 

The derivative of (4) with respect to A involves £^rj Sx = (£ • H)rjs x . The integral 
of this expression over S yields the first variation of the area. In terms of the mean 
value 




of (£ • H) over each leaf of the foliation, we can write d\S\\/d\ = a(X)\S\\. Using 
standard expressions for derivatives of geometric objects within an integral we now 
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find 



dM H {S x ) 



S 



A 



167T V 167T 



An X (S x )a(X) - 
+ 



Sx 



H 2 + -A 



(6) 



Note that only the normal part £ of v = d&(d\) appears in this expression because the 
tangential component of v gives rise to a divergence of a vector on T" which integrates 
to zero, S\ being closed. 

In order to reformulate (6), we recall some standard embedding formulas, e.g. [18], 
and collect some computations in Lemma 1 and Lemma 2 below. Let X, Y be tangent 
vector fields to S, £ be a normal vector to S, and V the covariant derivative on V. We 
can decompose the covariant derivatives along S into components tangent and normal 
to S as follows: 



DfY-K[X,Y 



(7) 
(8) 



Here D is the Levi-Civita connection on (S,h), K is the second fundamental tensor 
of S, V" 1 is the connection on the normal bundle of S and Ag : T p S — > T P S is the 
Weingarten map, both defined by the decomposition in (8). 

Lemma 1 Let S be a C 2 embedded submanifold of arbitrary codimension in a Rieman- 
nian manifold (V, g) of arbitrary signature and assume that S has a non- degenerate 
induced metric. Let X,Y be vector fields tangent to S and £ a C vector field normal 
to S. Then 



(£,K)(X,Y) 



± 



Proof. By construction d&\(X) is a vector field tangent to S\. It follows directly that 
the commutation of the velocity vector v = d$(d\) and X gives a vector field tangent 
to S and the same happens with [£, -X"], i.e. [^X]- 1 = 0. Using (7) and (8) and the 
definition R(a, (5)^ = V^V^ — V^VsJ — V[-^7 for the Riemann tensor we get 



+ k(y, [£,X]) = 
(V^f - V ( K(X, Y) - V^V f i- [I Y]Y + K (?, [£ X] 



+ 
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Using now Vg(K(X, Y)) = Vg { ^^+£^K(X, Y)) and the fact that £g(K(X, Y)) = 

£^(K)(X, Y) + K([l X],Y) + K(X, [£ Y]), the Lemma follows directly. □ 

We now treat the non-trivial term £^H in (6) as follows. As H is normal to S, 

only (£^H) ± needs to be calculated. We choose a basis {e^} of the tangent space Tj. 
Without loss of generality, we can assume that this basis is Lie propagated along the 
variation vector £, i.e. [£, e A ] = 0. Indices A,B,--- therefore refer to tensor objects 

— * — * 

within S expressed on this basis. For instance Kab = -^(eU, ee), and Va = V^* A . 
Such indices are lowered and raised with the induced metric Hab and its inverse h AB . 
Directly from the definition of 



H = tr Sx K = h AB K AB 



(9) 



and d x h AB 



-2{K AB • £ ), we get, after making use of Lemma 1, 



(£&) 



2k ab (k ab ■ e 

h AB K[A^eA),e B )-h AB 



R(i,e A )e B \ ± - (V^f)" 1 . 



(10) 



The first and fourth terms are the same apart from coefficient, while the second and 
third terms combine as — tr Sx ^V ± V- L |^. Also, inserting (10) in (6) we observe that 

the last term, in V^£, cancels with (£^g)(H, H). Notice that this term is not intrinsic 

for the variation because it depends on how £ changes along normal directions. It is 
clear that such dependence cannot occur in a first derivative. Summarizing, we have 
the following result. 

Lemma 2 With the same assumptions as in Lemma 1, suppose also that S is closed 
and orientable. Then, the variation of the Hawking mass along £ reads 



dM H (S x ) 
d\ 



1 



+ 



4vr X (S A )a(A) + 

l07T V l07T L 

2(K AB ■ H)(K AB ■ £) - 



AB 



(£■ h) + i (A) 



H 2 + -A) + 



+ 2h AB e A )e B ■ H + 2tr Sx H • V^V^f rj s . 



:n) 



We now restrict ourselves to spacelike two-dimensional surfaces in a four dimen- 
sional spacetime of Lorentzian signature. In particular, the normal spaces are now 
two-dimensional and Lorentzian, and inherit an orientation from the orientation of S. 
We can introduce projection operators P^ v and P^ u which act as the identities on 

Tp and T^, respectively, and annihilate vectors in the other orthogonal spaces. In 
particular, we have gw = + Pf. 

— * 

This setup also allows us to define the Hodge dual operation on normal fields W, 

1 



W* 



^sW^y B rj AB 
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We have (W**) = W and (W ■ U*) = — (W* ■ U). Therefore any normal vector is 
orthogonal to its dual and both have opposite norms. In particular, a normal null 

— * — * — * — * 

vector is either self-dual, I* = I, or anti-self dual, I* = —I. We note the following 
completeness property, valid for any pair of vectors u, v G T^: 



v a vP + uV - v* a u* p - u* a v^ = 2{v-u)Pf. 



(12) 



This expression can be proven directly by noting that both terms are symmetric, or- 
thogonal to S\ on each index and give an identity when contracted with any tensor 
product of u a , v a , w* , v *. In particular, when k and / are null and linearly independent, 

— * — * 

with = — (/ • k) 7^ 0, we have 

k a f + rk 13 = -<f>Pf. (13) 

— * — # 

We next introduce the Ricci scalar R{h\) of (S\, h\), and the trace-free part H — K — 

-. — * 

^Hh of the extrinsic curvature vectc 
part of the gravitational energy" by 



\Hh of the extrinsic curvature vector ^.Furthermore, we define 6^, the "transverse 



tAB 



--(U AB -U 



AB 



= (ti AB ■ h*) (? • n AB ) - \ (uab ■ n AB ) (f • h* ) (14) 



for any pair of vectors H, £. This terminology will be justified later, and the equivalence 

— * — # 

of the two alternative expressions in (14) follows from (12) with v — > H and u — > £. 
We also define 

n = -R(h x ) - -H 2 - -A (15) 
2 v x! 4 3 v ) 

which may be called the "Hawking energy density" since, using the Gauss-Bonnet 
formula, 

(16) 



we find from (4) that 



M H (S) 



[ R(h x )rj Sx = 4n X (S x ) 



is ,1 



167T \47T J s 



(17) 



Lemma 3 Let (V, g) be an oriented 4- dimensional spacetime of signature (— , +, +, +) 
and S a 2- dimensional oriented closed and spacelike C 2 surface in V . Let $ denote a 
flow of two-surfaces as defined above. Then, the derivative of the Hawking mass along 
the flow reads 



8ir V 167T J s 



dM H (S x ) 
dX 

+ tr Sx \H • V x V x ej - ft 

where G is the Einstein tensor of (V,g). 



(C(H\ f ) + A(H* • f )) + 8irO T (H*, t )+ 



(f H) - a(X) 
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Proof. We first note that 

R a pH a e = Ra(sH* a e P + (f - H) PfRcfi. (19) 

which follows from (12) with v — > if and u — > £ by contracting with i? Q /j. We next 
recall the Gauss identity 

iflfR^ps = R(h x ) -H 2 + (k ab ■ K AB ) . (20) 

We can now show that 

h AB (i?(£, e A )e B ■ H) = G{H\ £* ) + \ (f - &) [h 2 - R(h x ) - (k ab ■ K AB )] (21) 

by decomposing £*, if and P± u in a null basis k, /, using (13) for P^ and working out 
the contractions, using also (19) and (20). Finally, we obtain (18) by inserting (21) 
into (11) and making use of the Gauss-Bonnet formula (16). □ 
We now reformulate (18) in the case where £ is not null, and we write (£ • £ ) = £ 2 = 
ee 2i> f or e — ±1 anc i some function ip. We decompose £ in a null basis as follows: 

£*=£; + £jt where £/ = Al, £fc = Bk for some functions A, B. (22) 

Then ee 2 ^ = — 2AB<p ^ 0. Furthermore, we introduce the following one-form on S\ 

1 (i-v A £ k k-v A i 



u 



A 



20 I B A 



and the "longitudinal part of the gravitational energy density" L (£) MJ , 

8nQ L (£)(H*,e ) = (\U\ 2 + W| 2 )(£- H) - (2U ■ EhJ>)(£ ■ if*). (23) 

— * 

Lemma 4 We require the assumptions of Lemma 3 and in addition that £ is nowhere 
null. We obtain 



dM H (S x ) 1 \S. 



J [(g(H\ £* ) + A(if * ■ £* )) + 8vre T (ff *, £* ) + 



G?A 87T V 167T 7 S 

+87re £ (f)(H*,f ) - D.U{£- if*) + [(£ • if) - o(A)] (A^ - n)] rj Sx (24) 

where D.U denotes the divergence of U . 

Proof. In terms of the basis introduced above, we have 

V x £ = (Dil> + U% + (Di/> - U%. 

It follows that 

tr 5A (V ± V ± £) = (Aip + D ■ U + \Dip + U\ 2 ) & + (Aip - D - U + \Dip - U\ 2 ) £ fe 
and hence 

tr 5A (if • V ± V ± £) = (Ai/> + \Dij\ 2 + |f/| 2 ) (£ • if) - (D ■ U + 2U ■ D^) (£ • if*). (25) 
Inserting this in (18) and using (23) gives (24). □ 
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Causal Past 



Figure 1: Decomposition of the normal space into four closed sets. 



3 Monotonicity properties of the Hawking mass 

Lemmas 3 and 4 are the basis for a systematic discussion of monotonicity of the Hawk- 
ing mass. In both expressions there are four types of terms. We shall discuss them 
separately and show that each of them is non-negative under suitable conditions on 
the flow. Finally we put these conditions together to obtain the monotonicity result 
for Mh-i Theorem 1. 
As to the integral 



G(^,f)+A(^-f)]r/5 A (26) 



we require that 



{G aP + Ag afi )u a v^ > (27) 

for any pair of future (or past) causal vectors. If the Einstein field equations hold, the 
bracket is equal to 8-7T times the energy momentum tensor. In any case, we call (27) the 
dominant energy condition. Accordingly, (26) is non-negative if H* and are future 
(past) causal. 

In order to consider the following terms in (18), it is convenient to define four 
closed subsets of the normal space T^. Being a Lorentzian vector space, we can speak 
of causal (i.e. timelike or null) vectors and of achronal (i.e. spacelike or null) vectors. 
Let us denote by the subset of future (past) causal vectors and the set of vectors 
such that its Hodge dual is future (past) causal. Notice that all the vectors in At- are 
achronal. These four sets clearly cover the normal space and they are not disjoint. 
The vector belongs to all of them. Furthermore, the intersection C + fl A + , dropping 
the subindex p from now on for simplicity, is the set of self-dual future null vectors, 
and similarly for C~ fl A~. As for C + fl A~ , this is the set of anti-self dual future null 
vectors, and past for C~ fl A + . This decomposition is visualized in Figure 1. 
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A positivity Lemma for T (i/*,£*) follows directly from the positivity result for 
the so-called super-energy tensors, Theorem 4.1 in [19]. Since our case only requires a 
simpler version of the result, we add a proof for completeness. 

Lemma 5 The inequality 

2 (f AB ■ H) (f AB ■ i) - (f AB ■ f AB ) (£-H)>0 (28) 

holds for any tensor Tab if and only if£, H G C + , £,H G C~ , £,H G A + or £,H G A - . 

Proof: The "only if" part of the Lemma is easy to prove by finding suitable coun- 
terexamples. For the direct part, let p be any point in S. The object 2T AB ^T iy 



AB 



^AB^f B, y^\p defines a map on T p N x T P N which is obviously continuous. Thus, it 
suffices to prove the inequality almost everywhere on the subspace (C + x C + ) U (C~ x 
C~) U (A + x A + ) U (A - x A - ). We can assume without loss of generality that H 
and C, are both non-null and linearly independent. Since, by assumption, both vec- 
tors belong to the same subset C + , C~, A + or A~ , we have iJ 2 = ea 2 , £ 2 = e& 2 , 
where e = ±1 and a, 6 are strictly positive. Moreover we have e(£ • H) > 0. On any 

two-dimensional Lorentzian space, the inequality (w • v) 2 > v 2 u 2 holds for any pair 

— * — » — » — * 
of vectors. Hence e(£ • H) > oh. Since /J and £ are linearly independent, we can 

— * — * — # 

decompose Tab = ^abH + f^AB^- Direct substitution into the left hand side of (28) 
gives 

e (£■ H) (S AB S AB a 2 + Q AB n AB b 2 ) + 2a 2 6 2 S AB fi AB > ab (a 2 S AB S AB 
+b 2 Q AB Q AB + 2a6S AB fi AB ) = ab (a^ AB + 6ft ab) + bQ AB ) > 

and the Lemma is proven. □ 
Putting T AB = H AB , it follows that T (i7*,£*) is monotone iff the velocity vector 
of the flow points into the same quadrant as the mean curvature vector, at each point 
on the surface, or equivalently if £, H G C + , £, H G C - , £, H G A + or £, if G A - . 
We turn now to the positivity of Q L (^)(H*,^). This object is defined only when 

— # — # 

£ is non-null. However, equation (29) will be used later in the limit when £ tends to a 
null vector. 

— * — * — * 

Lemma 6 Let £ be non-null everywhere on S\ and asume that either (i) £ G C , H G 
C~ ; ^ £ G C~,H G C + ; (raj £, H G A + ; or (raj £,HeA~ holds everywhere on S\. 
Then Q L {£) (H*,£*) > 0. 

— * 

Proof: Conditions (i), (ii), (iii) or (iv) together with the fact that £ is not null 

— * — * 

implies (£ • H) ^ everywhere. We rewrite (23) as 



87r0 L (O(^,f ) = (i-H 



2 



Tj2.pl 

+ ^V|^| 2 (29) 



after using H 2 £ 2 = (£ • H) 2 - (£ • iP) 2 , which follows from (12). Each of the hypotheses 
(i) to (iv) implies H 2 ^ 2 > and (H ■ £) > 0, and the Lemma follows. □. 
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Turning now to the other terms in the derivative of the Hawking mass (18), we 
first assume that that £ is non-null everywhere, in which case we can use (4). The last 
term, namely 



involves the product of two factors. We first note that the geometric flow of 2-surfaces 
is obviously independent of the parametrization used, and a reparametrization changes 
the velocity field £ by a nowhere-zero constant on each leaf. In particular, such a 
reparametrization rescales a(A) by an arbitrary nowhere-zero function of A. Thus, if 
a(A) 7^ everywhere, a = 1 can be assumed without loss of generality. If, on the other 
hand, a has zeros, be it at isolated points or on intervals, those values are meaningful 
and independent of reparametrization. 

— * — * 

The first factor in (30) is the difference between (£ • H ) and its mean value on 
S\, a(A). Consequently, it will never have a constant sign on the surface unless it is 
identically zero. While it might be possible to define a flow in such a way that each 
factor changes sign at the same place, such a condition would be very complicated to 
analyze. Hence we restrict ourselves to the cases that either the first factor is identically 
zero or that the second one is constant, each of which remove the integral (30). 

As to the first option, the condition (£, H) = a(A) is the generalization to codimen- 
sion two of the IMCF on surfaces in a three-dimensional space. Indeed, in codimension 
one, the mean curvature vector and the normal velocity are obviously parallel, and 
setting a — 1, assuming a(A) ^ 0, this condition just states £ = fl/p, where fx is the 
unit normal vector and p the mean curvature, as before. This is precisely the IMCF 
condition. Notice, however, that IMCF in codimension one fixes uniquely the velocity 
of the flow while in codimension two, it still leaves room for imposing one extra scalar 
condition on the velocity. We shall use the name IMCF in both cases. 

The requirement which could be used alternatively to remove (30) reads Aip — Q = 
a(X) for some constant a(X). By Gauss' law and by a trivial Fredholm argument, this 
equation has a unique solution for ij) iff a(X) = j x Qr)s x - Curiously, this integral is 
related to the Hawking mass itself, c.f. (17). As for the IMCF condition, this one has a 
counterpart in codimension one as well, obtained by substituting the mean curvature p 2 
for H 2 in f2, c.f. (15). This condition could be used in place of the IMCF to guarantee 
monotonicity of the Geroch mass. The question whether the corresponding flow exists 
will be discussed elsewhere. 

In the case where the £ is allowed to be also null, the splitting (24) is not possible 
and we must return to (18). The last term in this expression can be made zero with 
sufficient generality only by imposing the IMCF condition (£ • H) = a(A). 

It only remains to consider the integral 




(30) 




(31) 



Lemma 7 Under the same hypotheses of Lemma 3. 

— * — * —. 

(a) If H\s(\) = B£*\s(\) for some function B on the surface, then 1 = provided £ 
is null all over S\ or if B is constant on S\. 
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(b) Assume (£ ■ H) = a(X) and (£ • if*) = c(A) for some constant c(A) on S x . Then 
I > provided (i) £eC + ,He C~ , (ii) £ G C~,H G C+, fmj £ ', H G A+ , or 
(iv) £, ff G A - /ioWs everywhere on S x 

Proof: Regarding part (a), if is null on an open set U C S x , we in fact have if = ±-B£* 
on U and hence Vi^ = &a£ on [/, for some b&- Thus, V^V^^is parallel to and 
hence orthogonal to ff. Consequently (H ■ V V £)\u — holds irrespective of the 

— * — * — * 

values of B. If £ is non-null everywhere then from the fact that (£ • if ) = and 
(£ ■ if*) = eBe 2lp , (25) implies that the integrand in i can be written as -BD A (£ 2 U A ) 

— * 

which integrates to zero provided B is constant. If £ is non-null almost everywhere the 

— * 

same conclusion holds by continuity. Moreover, if £ is null on an open set, the integral 
on that set is zero for any B, as shown before. Thus if B is constant on S x the integral 
vanishes irrespective of the causal character of £. 

Turning to (b), we first notice that if £ is non-null everywhere, then (25) to- 

— * — * — * — * 

gether with the fact that (£ • H) and (£ ■ H*) are constants on S\ implies I = 
J s ^ <d L (£)(H*,£* )rj Sx > 0, where Lemma 6 has been used for the last inequality. To 

— # 

include the case when £ may be null we use a continuity argument. The integral I 
defines a map from C°(S\) x C 2 (S\) — > R, where the first factor refers to H and 
the second factor to £. This map is obviously continuous with respect to the supre- 
mum norm on these spaces. Thus, it is sufficient to observe that the inequality holds 
for H,i e mt[C+(S x ) x C 2 (S X )] U intfC^) x C+(S X )] U int[A+(S A ) x At(S x )} U 
int[AQ (S\) x A 2 {S\)\ where Cq(S\) denotes the set of vectors h G C°(S\) such that, 
for all p G S\, h(p) G C+. C^Sx) is defined analogously by replacing C°(S\) by 
C 2 (S\) and the other sets are defined similarly. Notice that this open set is dense in 
[C+(S X ) x C-(S X )} U [C-(S X ) x C + (S X )} U [A+(S X ) x A+(S X )] U [A~(S X ) x A~(S X )], 
where we want to prove the Lemma. □ 

— * — # 

Remark 1. Condition (a) holds for instance when (H ■£) = everwhere on S x and 

— # 

£ does not vanish on open sets of S x . 

Remark 2. This Lemma tells us that in the generic case, i.e. (£ -H) not identically 
zero, we need to restrict the scalar product (£ • H*) to being constant on each leaf. We 
call this the dual inverse mean curvature condition. Such a condition is unnecessary 
in the time-symmetric case, where the Geroch-Hawking mass depends only on the 
intrinsic curvature of the hypersurface generated by the flow. 

Among the conditions (i) to (iv) that ensure monotonicity in part (b) of Lemma 
7, only (iii) and (iv), i.e. velocity vector and mean curvature vector both achronal 
and to the same quadrant, give conditions which are compatible with those ensuring 
non-negativity of (26) and G T (if*, ). 

— * — * 

Regarding the case a(A) = on S x , the orthogonality of H and £ implies that they 
cannot belong to the same quadrant, thus violating the non- negativity condition in 

— * — * — * 

Lemma 5, unless £ is null and H = B£, with a non-negative arbitrary function on S x . 
Thus, when the surface is marginally trapped and the flow is null, monotonicity is very 
easily achieved initially, the only condition being that the velocity is pointwise parallel 
to the mean curvature of the surface. Of course, the property of H being null will not 
be maintained by the flow and the flow vector will have to adjust, if at all possible, to 
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keep a monotone flow. 

Summarizing all the results above and noticing finally that the integral f s D.U{^ ■ 

— * 

H*)t] Sx in (24) vanishes provided the dual IMCF condition is satisfied, we can write 
down the following theorem which gives sufficient conditions on a spacetime flow of 
2-surfaces in a spacetime so that the Hawking mass is monotone along the flow. 

Theorem 1 Let S be an oriented, spacelike, closed, C 2 embedded 2-surface in a space- 
time iy,g). Let £ be a C 2 normal vector field on S, H the mean curvature vec- 
tor of the surface and M H (S) its Hawking mass. Let S\ be any flow of spacelike 
two-surfaces starting at S with £ as normal component of the initial velocity. Then, 
dM H (S\)/d\\\ =0 > holds whenever the following four conditions hold 

(1) The spacetime satisfies the dominant energy condition (27). 

(2) The mean curvature vector H is spacelike or null on S, and £ points into the 

— * 

same causal quadrant as H . 

(3) Either (a) the IMCF condition holds, i.e. (£ ■ H) = a®, where a$ is a non-negative 
constant, or (b) £ 2 ^ everywhere on S and 

= ±R(h x ) -±H 2 + h-a (32) 

where a is constant. 

(4) The dual IMCF condition holds, i.e. (£ • H*) = c , where c is a constant. 

Remark. We now assume that £ 2 ^ and that (3), (4) hold. Eqn. (24) can be 
reformulated so that the energy interpretation becomes transparent. We can write the 
change of Mjj as follows 

dM H (S x ) f ( 1 



d\ 



Here we have introduced x — H* ^\S\\/lGir; the normalization is chosen such that 
X coincides with the asymptotically unit timelike Killing vector when S is a 2-sphere 
on a time-symmetric slice in Schwarzschild. In more general situations with timelike 
Killing vectors, we are not aware if they necessarily coincide with x- We note that on 
any marginally trapped surface, if* and therefore x is a nv ^ vector. A spin-coefficient 
version of the energy conservation law (33) is given in [20]. While Ot is a tensorial 
object (14) depending only on the geometry, <d L depends on the velocity vector £ as 

— * — * 

well. In any case, the components of these objects with respect to a null basis (k, I) 
with (k • I) — —(p read 

SnOl = (l-U AB ) 2 8nO T kk = (k ■ U AB ) 2 B T H = Of k = 0. 
e L (0« = Q L (£) kk = %TtQ L {£) lk = <t>{Di)-U) 2 MQ\i) kl = <PW + U) 2 . (34) 
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These tensors may be interpreted respectively as encoding transverse and longitudinal 
modes of the gravitational field, with 6 = 6 T + L taking the same form as in a 
corresponding energy conservation law along black-hole horizons [21]. In particular, the 
components and G^, can be interpreted as energy densities of ingoing and outgoing 
gravitational radiation, taking the same form as standard definitions at null infinity in 
an asymptotically flat spacetime. In that case, the conservation law corresponding to 
(33) is the Bondi mass equation [22]. 

We proceed by discussing the requirements of Theorem 1 in some detail, restricting 
ourselves to the IMCF condition, i.e. (3. a). 

First a remark on terminology: recalling that £^rjs x = (£ • H)rjs x , we see that the 

— # 

first of these conditions implies that the area element is preserved under the flow £, 
while the second one implies the same along the dual £*. This motivates the terminology 
uniformly expanding flow (UEF) [6] for both (£ • H) — a(X) and (£ ■ H*) = c(A), while 
we reserve the term inverse mean curvature flow (IMCF) for the first one alone. The 
UEFs have some interesting special cases which we discuss in turn. 

We first assume that the flow velocity is null everywhere, while the mean curvature 
vector of the surface is achronal everywhere. To obtain monotonicity, we impose the 
IMCF condition with a(A) = 1. This implies the dual IMCF condition since = ±£, 
and therefore 

tf*) = - (f -h)=T (f H) = =Fl. 
All the conditions of the theorem hold provided £ and H point into the same quadrant. 

— * — * 

Notice that the IMCF condition actually prevents £ or H from vanishing anywhere on 
the surface. This null monotone flow was discovered by Hayward [6]. Since £ is null 

— * — * 

and hypersurface orthogonal, the flow is geodesic, i.e. V|-£ = F£ for some function Y. 
In terms of an afnnely parametrized geodesic with tangent /, we have £ = QzH and the 

1 ~~ * 

flow can be determined by solving the ODE £(A) = 9Z /(A) = 1. We emphasize that 
this is the only case for which we obtain local existence. 

We next assume that H is spacelike everywhere, i.e. H 2 > 0, which allows us to 
write the velocity vector in the form 

£ = H~ 2 [H - c(\)H*) , (35) 

choosing a(A) = 1. A special case is that c(A) vanishes identically, so that the velocity 
vector of the flow is £ = H~ 2 H. This flow was mentioned by Huisken and Ilmanen 
in [9] and analyzed in more detail by Frauendiener [17], who showed monotonicity in 
an explicit way. In the general case covered by (35), we can set a(A) = 1 so that the 
velocity flow contains one arbitrary function of one variable c(A) satisfying |c(A)| < 1. 
In the particular case |c(A)| < 1, this is the spacetime version of the flows found by 
Malec, Mars and Simon [16] in terms of initial data sets. Indeed, in this case the vector 
£ in (35) is spacelike by construction, so that the hypersurface £ generated by the flow 
of 2-surfaces will be spacelike, if it exists. Thus, when the flow is viewed as a flow of 2- 

— * — * 

surfaces in this 3-dimensional Riemannian manifold, the IMCF condition (£-H) = a(A) 
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translates into the condition that S flows within E as an inverse mean curvature flow. 
The other condition (£-H*) = c(A) becomes a condition involving p, the mean curvature 
of S C E, and g, the trace on 5 of the extrinsic curvature of E as a hypersurface in 
spacetime, which explicitly reads q/p = c(A). This is precisely the condition obtained 
in [16] (see also [23]) in order to ensure monotonicity of the Hawking mass in an 
initial data setting. These flows significantly generalize flows along the inverse mean 
curvature vector, since the one-parameter freedom after fixing a(A) = 1 allows flows 
from a given initial surface to cover a spacetime region. In particular, while H/\H\ 2 

— * 

can only approach spacelike infinity in an asymptotically flat spacetime, £ of the form 
(35) can also approach null infinity. This suggests the possibility to prove the Penrose 
inequality (1) not only for the ADM mass M but also for the Bondi mass Mb < M. 
The latter inequality is significantly stronger than the former, but it also follows from 
the heuristic argument involving cosmic censorship [1]. 

The final special case of the UEF guaranteeing monotonicity of the Hawking mass 
are surfaces with null H and velocity vector which is either spacelike or null, towards 
the same quadrant as H. Since H = ±H*, the dual IMCF condition implies the 
IMCF condition, as in the case where £ is null. The surface S is restricted to being 
a future (past) marginally trapped surface, which is precisely the starting surface for 

— * — * 

flows aiming at proving the Penrose inequality. Given H, the velocity vector £ at S has 
a freedom of one arbitrary function. Thus, we are allowed to flow out of a marginally 
trapped surface along a large class of flows for which Mh is monotone. 



4 On existence of uniformly expanding flows 

In this section we discuss the prospects of obtaining existence of the flows defined 
above. Recall that the method of Huisken and Ilmanen in codimension 1 rests mainly 
on two properties of the IMCF: firstly, the embedding functions determining the flow 
satisfy a parabolic system. Secondly, the level sets of the flow satisfy a degenerate 
elliptic system. As to codimension 2, Huisken and Ilmanen noticed that the embedding 
functions describing a flow along the inverse mean curvature vector, i.e. (35) with c = 0, 
satisfy a backward-forward parabolic system. There is no general theory which would 
guarantee local existence of solutions to such systems, and counter-examples indicate 
that such systems have to be treated with care. Unfortunately, as we will see below, 
the system keeps this vicious behaviour for any choice of c. 

Nevertheless, the more important component on the work of Huisken and Ilmanen 
is the weak (variational) formulation of the level set formulation for the flow. We have 
obtained a generalization of this formulation which we will describe in turn. 

We first give the description in terms of embedding functions $ Q (A, x A ) for a 2- 
surface S embedded in spacetime, following the notation introduced in Sect. 2. In 
local coordinates we have £ M = <9$ M /<9A and the tangent vectors to S can be written 
as e A M = d^/dx A . For the mean curvature vector we obtain 

r)(b v r)(bf 

AB ^ 11 1 up\ x=9 a q^a Q£B' y 6K) > 
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where r M are the Christoffel symbols of the ambient space, and A is the Laplacian 
on S. In this way, (35) becomes a system of PDEs of the form 



To determine the type of this system we need to calculate the eigenvalues of the sym- 
metrized part Q + of the matrix 

QTb = ~ z A z B (38) 



y dx A dx 

where z A is a unit vector in tJ, and f denotes the transpose with respect to some 
positive definite metric. Restricting ourselves now to codimension two, we obtain for 
the linearization of the operators H and H* acting on $, 

, r z A z B = -P ± % — — v z A z B = -rf s \ (39) 

° \dx*dx» ) ° [dx^dx^ J 

where we have employed the dualized 2-form 77^: 

1 „. 7 g A B 



It follows that 



Vs — 2^ /3~fS e A e B 7 ls 

2H"H P \ c(A) 



H 2 



(40) 



To do the symmetrization, we can choose the metric L a p = H 2 (H a Hp + H*H^) so 
that the symmetrized part of is independent of c(A) and reads 

(Q + QT P = ^ + H*°H$) (41) 

This matrix has H a and fP a as eigenvectors, with eigenvalues 2/H 2 and 2/H* 2 = 
—2/H 2 , respectively. Therefore (Q + Q^) a ^ is indefinite and (35) is a so-called forward- 
backward parabolic system. The previous reasoning applies in codimension one if all 
H* are removed, and yields parabolicity, as well known. It can be easily checked that 
the reasoning above is independent of the metric L^ v used for symmetrization. 

Turning now to the level set formulation, we start with recalling the case of codi- 
mension 1. One can envision the family of surfaces 5(A) as level sets of a level set 
function u(x) defined on a spacelike hypersurface £ of the spacetime. If u is twice 
differentiable, the inverse mean curvature flow condition can be written in terms of u 
as a degenerate elliptic equation, namely 

V. ( ) = M (42) 
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where V is the covariant derivative on E, and V.W denotes the divergence of some 
vector W. The variational formulation reduces the requirements on differentiability. 
While (42) are not Euler-Lagrange equations of a functional of u as the only variable, 
they can nevertheless be obtained by minimizing the functional 



with respect to variations of v on any compact set K. One can imagine this variational 
principle as a two-step procedure in which the term \T>u\ on the right hand side of (42) 
is first fixed or "frozen" while E u is varied with respect to v, and in the second step 
\T>u\ is put equal to the \T>v\ with the function v obtained from the minimization, with 
derivatives understood in a weak sense. 

As to codimension 2, it would be desirable to have some generalization of (42) with 
an extension of u to the extra dimension. While from the equation itself it is not clear 
how this could be accomplished, we turn directly to the variational formulation, which 
does, in fact, involve a natural extension of u. 

We start from some initial two-surface S embedded in some three-surface E. The 
latter is arbitrary at this stage, but will be fixed automatically by the variation prin- 
ciple. We consider the level sets S(X) C E given by the level set function u(X) of 
the IMCF near S, obtained via the Huisken-Ilmanen procedure. We denote by z/(A) 
and //(A) = ^*(A) the unit normals to 5(A) in the future timelike direction normal to 
E and in the spacelike direction tangent to E, respectively, assuming their existence. 
The variational principle should now move the flow in the direction of £ given by (35), 

— * 

i.e. E should have £ as tangent vector whenever the latter is defined. Equivalently, 
the mean curvature vector to S\ should be orthogonal to v + c(u)fi, and this latter 
property is used in the weak formulation. We extend the function u off E in some 
neighbourhood of E in such a way that the directional derivative of u in the direction 
v + c(u)/i, whose existence we may assume, is zero everywhere on S\. We denote by 
T any other spacelike hypersurface in this neighborhood which is equal to E outside 
some compact region T, and we take v to be a real-valued function on T which agrees 
with u outside of Y as well. We then define the functional 



In analogy with (43), the "frozen" structure now consists of the hypersurface E together 
with the level sets u, extended as above to a neighbourhood of E. The objects which 
are varied are T and v, and in the second step, E and u are adjusted according to the 
results of the variation. This justifies the following. 

Definition 1 Given c(A) G [—1, 1], then (E, u) is a weak solution to the corresponding 
uniformly expanding flow in an open region Q if it is a critical point of the functional 
Ej< (•, •) for all compact T C Q. 

We will also say that the 2-surfaces u are weak solutions to equation (35). Since 
u(x) may have constant regions, this allows the family of surfaces defined by the level 
sets of u(x) to jump, just as in the Huisken-Ilmanen case. 




(43) 




(44) 
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Given this weak definition, we may check that the above weak solutions agree with 
the usual smooth solutions when they exist. We first use the fact that (E, u) is a critical 
point of E£ uc (-, •) with respect to variations of E. Consider a variation T of E which 
is a compactly supported bump function on E times v + c(u)fi. Note that u does not 
change to first order in these directions. Hence, to first order, it is still true that 

E = / (\Vu\+u\Vu\) , 
irnr 

where we may extend u to be constant in this variation since we are at a critical point. 
But by the co-area formula, if we choose Y to be the region where a < u(x) < b, 

E= ( {l + \)\S{\)\d\, 

J a 

where (S^A)! denotes that area of 5(A). Since we are at a critical point, the first 
variation of E must be zero, which, in the case that 5(A) is a smooth family of surfaces, 
follows if and only if our variation direction is orthogonal to the mean curvature vector 
of S(\). Equivalently, it follows that £ c from equation (35) must be tangent to E. 

What remains to be done is to check that the actual flow vector, which is also 
tangent to E and orthogonal to each S(X), has the same length as £ c . Again we use the 
fact that (E, u) is a critical point of (■, •) with respect to variations of u(x). This 
implies that u(x) is a weak solution to inverse mean curvature flow in E in the sense 
defined by Huisken and Ilmanen. Thus, the level sets 5(A) of u(x) have the uniformly 
expanding property that the rate of change of the area form of the surfaces equals the 
area form itself. This proves that the actual flow vector of the surfaces is not only in 
the same direction as £ c but also has the same length, and therefore they are equal. 
Thus, the smooth flow case agrees with this weak definition. 

To conclude this exposition we wish to emphasize that having a definition of a 
weak solution to equation (35) is still very far from an existence theorem. However, 
the correct notion of a weak solution is a prerequisite to an existence theorem along 
the lines of Huisken and Ilmanen. 
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